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1, Introduction 

There are several studies on the simultaneous determination of cloud optical 
thickness and mean droplet size as derived from multispectral radiometers 
(Twomey and Cocks, 1982; Curran and Wu, 1982; Arking and Childs, 1985). 
These studies are very important to elucidate the cloud microphysical struc- 
ture. As pointed out by many investigators (e. g. Rossow et al., 1985; Wu, 
1985), however, many factors bother the successful retrieval when we apply 
these methods to real data. Therefore, many ground truth experiments are 
required in order to test the validity of these methods. Here, we have ana- 
lyzed the data taken in the marine stratocumulus intensive field observation 
component of the First ISCCP Regional Experiment (FIRE), conducted off the 
California coast in July, 1987 using the Multispectral Cloud Radiometer 
(MCR). 

The MCR has six near-infrared channels and one 10.7 pm-channel (Curran et 
al., 1981). The first channel (X = 0.754 pm) has proven to be very useful for 
estimating the cloud optical thickness (King, 1987). AddiAg the additional 
two wavelengths ( X = 1.645, 2.160 pm), we are able to estimate the mean 
droplet size (Curran and Wu, 1982). Arking and Childs (1985) also studied a 
similar problem using the 0.73 and 3.7 pm channels of the NOAA/ AVHRR. 
In spite of these pioneering works, there seem to be insufficient systematic 
studies of the efficiency of the algorithms. For this purpose, we have per- 
formed numerous simulations of the reflected radiation from clouds using 
an improved radiation calculation algorithm (Nakajima and Tanaka, 1988). 
Results of the simulation and an application of our algorithm to the MCR 
data are presented in Sections 2 and 3, respectively. 

2. Retrieval algorithm 

For the purpose of analyzing the three MCR near-infrared channels outside 
water vapor and oxygen absorption bands (X= 0.754, 1.645 and 2.160 pm), we 
investigated the characteristics of the numerically simulated bidirectional re- 
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flectivity of plane-parallel earth-cloud systems. The variables" we examined 
were the optical thickness (x) and the mode radius of cloud droplets (r m ), de- 
fined by the following log-normal cross-section spectrum: 

= C Q expj-(lnr — lnr m ) 2 / (2a 2 )j, (1) 

where C is the geometric cross section of droplets; r is the particle radius; o is 
the log-dispersion of the spectrum. Investigating the first and second mo- 
ments of several model size spectra of cloud droplets, it is found that a can be 
fixed at 0.35 for most problems even including aerosol polydispersions. The 
following points are clarified as a results of the numerical simulation: 

(1) As a first approximation, the cloud optical thickness x at X- 0.754 pm can 
be determined from the reflected radiation at this wavelength. The re- 
trieval is stable, because the cloud reflectivity increases monotonically 
with increasing optical thickness and the dependence on the droplet size 
is relatively small. 

(2) Using the cloud optical thickness determined in step (1), the cloud droplet 
mode radius r m is estimated from an absorbing channel, i.e., 1.645 or 2.160 
pm radiation. The optical thickness dependence is small as absorption 
increases. 

(3) Iterating steps (1) and (2) or investigating these steps for several pre-as- 
sumed particle radii, we can have an optimum set of x and r m . We have 
developed two methods for the simultaneous retrieval of cloud optical 
thickness and mode radius, a bispectral method and a multispectral 
method. The bispectral method uses the two channels of 0.754 and 2.160 
pm. This method is hardly affected by the vertical inhomogeneity of the 
cloud microphysics, but tends to be marred by noise in the absorbing 
channel. Also, this method can have multiple-solutions of r m , as 
demonstrated in Fig. 1. For simplicity we show in Fig. 1 a mode radius 
dependence of the cloud spherical albedo at X = 2.16 pm for x = 8, 16 and 
32. Figure 1 shows that the profiles are a peaked function of mode radius 
with a peak around r m = 2-3 pm. Thus, for example, a measured re- 
flectance at 2.16 pm, given an optical thickness derived from the 0.754 pm 
channel, could be produced by one of two possible mode radii. Therefore, 
it is very important to set a suitable minimum droplet radius for the 
algorithm. On the other hand, the multispectral method uses all three 
channels. The optimum value of r m is determined in this method so as to 
minimize the following root mean square deviation between the ob- 
served and theoretical values: 




where n = 2 stands for the second channel (1.645 pm) and n = 3 for the 
third channel (2.160 pm); superscripts c and m stand for theoretical 
(computed) and measured values, respectively. In the above equation, x 
is determined from the first channel (0.754 pm), as in the bispectral 
method. This method is relatively insensitive to noise in any individual 
channel and we can have a unique solution for r m , potentially different 
from the bispectral method. But this method fails in some case of vertical 
inhomogeneity. We have adopted both methods for the sake of diagnos- 
ing the soundness of the retrieval. 

(4) The error involved in the retrieval increases in the backward portion of 
the scattering. The error in r m increases with decreasing r m ; the errors in 
x and r m increase in the region x < 6 or x > 100. 
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Fig. 1. Spherical albedo of a plane-parallel cloud versus the droplet mode radius at X= 2.16 
pm. 

3. Analysis of MCR-data 

Since the MCR can produce a 2-D image of the upwelling radiation, we tried 
to develop a retrieval scheme computationally very efficient. We have 
generated four tables of bidirectional reflectivity for the scaled cloud optical 
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thickness x* « (1 -g)x = 0.4, 0.8, 1.2 and «>, where g is the cloud asymmetry fac- 
tor. Using asymptotic theory for x* 5: 1.8 (King, 1987) and computations for x* 
< 1.8, we can interpolate the bidirectional reflectivity with errors less than 2%. 
Since our intention is to derive effective parameters to describe the radiative 
characteristics of clouds, we assume homogeneous and plane-parallel stratifi- 
cation for the analysis. According to numerical simulations, x can be esti- 
mated with an error less than 20% and r m can be estimated at the effective 
depth of radiation unless the optical thickness and particle radius are too 
small. 

Figure 2 shows contours of the cloud optical thickness derived for the first 
half of the fourth ER-2-flight line of July 10, 1987 (about 1000 scan lines, 60 
km), A photographic print of the original data for this example shows that a 
band of x-15 was the edge of a large flat cloud deck spreading ahead of the 
flight course. The optical thickness varied from 5 to 20 around this edge. 

Figure 3 shows the horizontal distributions of x and r m along a constant pixel 
line (solid lines for pixel number 84 and dotted lines for pixel number 252, 
where the total number of profiles in an active scan is typically 336). 
Although r m showed a small monotonic decrease from left to right, the value 
was remarkably constant with values around 15 pm, in spite of the large 
change of x from 5 to 15 around the horizontal distance of 30 km (solid line) 
and large hole around 40-50 km (dotted line). Although we showed the 
value of r m derived by the multispectral method, there was no large difference 
from the bispectfal method except for small optical thickness < 6. For small 
optical thickness, however, there were some cases when the algorithm 
selected the smaller branch of the solution depending on the viewing 
direction, as shown by several steep dips of the mean radius. Since we had a 
large difference between the results of the multispectral and bispectral 
methods in this case, there may be large uncertainties in the retrieved value 
of r m around these dips (as expected by item (4) in the preceding section). We 
show the minimized x 2 -value in the multispectral method in the uppermost 
section of Fig.3. The rms-difference between the measured and theoretical 
values was less than 5%. 

4. Concluding remarks 

We have installed an efficient retrieval scheme for deriving the cloud optical 
thickness and droplet mode radius; and we have described the behavior of the 
retrieval error. Generally the scheme can retrieve the optical thickness and 
mode radius adequately unless they are too small; the use of the forward scat- 
tering region is more sound than the use of the backscattering portion. 




Fig. 2. Cloud optical thickness distribution for the fourth flight line of July 10, 1987. 
Horizontal and vertical axes are for scan-line and pixel numbers. Bold solid and dotted lines 
shows the sampling lines for Fig. 3. 



Fig. 3. Cloud optical thickness (x), droplet mode radius retrieved by the multispcctral 
method (r m ), and x along the bold solid and dotted lines in Fig. 2. 




Applying our method to real data, we have derived the two-dimensional dis- 
tribution of optical thickness and mode radius for a portion of one of the FIRE 
marine stratocumulus missions. In this case study, the droplet size showed a 
more uniform distribution than optical thickness with some correlation be- 
tween large droplet size and small optical thickness. Although we can find 
microphysical reasons for these tendencies, we are suspecting that the re- 
motely sensed droplet size may be overestimated. As a future problem, we 
will compare our results with in situ data of the droplet size distribution. 
Also it will be very important to check several reasons why the droplet radius 
might be overestimated, e.g., soot contamination, effect of escaping photons 
from the lateral sides of broken clouds, and so on. 
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